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Abstract: We investigate n-gluon scattering amplitudes in the multi-Regge region of N = 4 
supersymmetric Yang-Mills theory at strong coupling. Through a careful analysis of the ther- 
modynamic bubble ansatz (TBA) for surfaces in AdS§ with n-g(lu)on boundary conditions we 
demonstrate that the multi-Regge limit probes the large volume regime of the TBA. In reaching 
the multi-Regge regime we encounter wall-crossing in the TBA for all n > 6. Our results imply 
that there exists an auxiliary system of algebraic Bethe ansatz equations which encode valuable 
information on the analytical structure of amplitudes at strong coupling. 
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1 Introduction 

The computation of gluon scattering amplitudes in gauge theories such as Quantum Chromody- 
namics (QCD) or its super symmetric cousins is a daunting task. Over the last few years much 
progress has been made in the context of N = 4 supersymmetric Yang- Mills (SYM) theory, both 
at weak and strong coupling. These exciting developments exploit new hidden symmetries, such 
as dual conformal symmetry [1], and the celebrated duality with string theory on AdS§ x 5 5 
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[2]. There is some hope to find expressions for the amplitudes that are valid for all values of 
the t'Hooft coupling, at least in the multi-color limit. 

Such hopes were first nurtured by the intriguing BDS formula of Bern, Dixon and Smirnov [3]. 
It encapsulates the known infrared and collinear behavior of n-particle maximally helicity vio- 
lating (MHV) amplitudes in the planar approximation. The authors of [3] conjectured the BDS 
formula to determine the amplitudes at each loop order L > 2, possibly up to some additive 
finite function of the kinematic variables, the so-called remainder function. Initially, 
was suspected to vanish, i.e. the BDS formula was believed to be exact. 

Both gauge and string theory arguments subsequently confirmed this suspicion for n = 4, 5. 
In the weakly coupled theory, perturbative computations uncovered the before mentioned dual 
conformal symmetry of scattering amplitudes [1]. It implies that the remainder functions R^ 
can only depend on conformal cross ratios, i.e. on conformally invariant combinations of the 
usual kinematic variables. Since there are no such cross ratios for n = 4,5, the corresponding 
remainder functions have to be trivial. In other words, dual conformal invariance predicts that 
the BDS formula is exact for n = 4, 5 to all loop orders. This prediction was confirmed by a 
string theory computation of the leading term at strong coupling [4]. We shall say a bit more 
about the string theoretic analysis below. 

On the other hand, the remainder function R^ is now known to be non-zero for n > 5 and 
beyond one loop [5, 6]. Several authors have described tests of the BDS formula that exclude 
a vanishing remainder function. One of the most direct ways to see that R^ ^ is based on 
a study of the SYM scattering amplitudes in the leading logarithmic approximation, see [7, 8]. 
The high energy (Regge) limit probes the remainder function near special points in the space 
of kinematic variables. While the Regge limit of the function R^ vanishes at some of these 
points, for example when the limit is taken with all energies negative, the authors of [7, 8] were 
able to identify one region in which the Regge limit of R^ is non-zero. Hence, R^ must be 
a non-vanishing function of the kinematic variables. The analysis shows how computations in 
the Regge limit can provide strong and highly efficient constraints on the remainder function 
and its analytical structure. 

In the meantime, the analytic expression [9] for the exact two-loop calculation of the six-point 
function [10, 11] was used to perform the relevant analytic continuation into the region with 
non- vanishing Regge limit [12]. The results are in full agreement with [8]. This settles the 
remainder function in the two-loop approximation, and it supports the all-order leading log 
generalization in [7, 8]. More recently, progress has been made with the extension of the calcu- 
lation of R^ to n > 2: in [13] the symbol of R^ has been determined (up to two parameters), 
and in [14] this result has been confirmed, fixing also the two previously unknown constants. 
In [15] the authors quote results for the symbols of Rq for four loops (again up to a number 
of unknown constants) [16]. The form of the scattering amplitudes in the Regge limit which, at 
weak coupling, was derived in the leading logarithmic approximation is quite general, and is ex- 
pected to be valid also outside the weak coupling limit. As a function of the energy variables the 
amplitude contains Regge cut terms with power like dependence on s-like kinematic invariants 
(see below). The exponents depend on the kinematical region and are determined by the lowest 
eigenvalue of the BFKL color-octet Hamiltonian for an n-gluon system. These eigenvalues have 
recently been calculated in NLO accuracy in [17], and in [15] in next-to-next-to-leading order 
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(NNLO). The power- like energy dependence of the scattering amplitudes is multiplied by Regge 
impact factors which are now known also in NLO [18] and even in N 3 LO accuracy [15]. A first 
generalization of the leading logarithmic analysis to the 7-point amplitude has been started in 
[19]. 

The BFKL color-octet Hamiltonian possesses the very interesting property that it coincides 
with the Hamiltonian of an integrable open spin chain [20] in leading order. Hence, the weakly 
coupled theory provides direct evidence for integrability in the high energy behaviour of planar 
scattering amplitudes. 

Having reviewed all these results from gauge theory it is natural to ask what string theory has 
to say about the high energy limit of the remainder function . In order to understand how 
the issue can be addressed, we need to briefly sketch the development that was initiated by the 
work [4] of Alday and Maldacena. The main insight of this paper was the identification of the 
leading contribution to an n-gluon amplitude at strong coupling with the area A n of some 2- 
dimensional surface S n inside AdS§ . According to the prescription of [4] , S n ends on a piecewise 
light-like polygon on the boundary of AdS§. The light-like segments of this polygon are given 
by the momenta pj of the external gluons. For n = 4 it is possible to find the surface explicitly 
and the resulting amplitude matches the prediction of the BDS formula. Constructing S n for 
n > 5, however, turned out to be a rather difficult problem, at least for finite n and generic 
choice of the external momenta. The issue was resolved through a series of papers [21-23] in 
which the area of S n is related to the free energy of some auxiliary quantum integrable system. 
More precisely, it was argued that A n may be computed from a family of functions Y a;S with 
s = 1, 2, 3 and a = 1, . . . , n — 5. The latter can be determined by solving a set of coupled 
non-linear integral equations. Very similar mathematical structures are familiar from the study 
of ground states in 1-dimensional quantum integrable systems on a circle of finite radius R. 
Moreover, the functional A n resembles expressions for the free energy of such systems. So, in 
the sense we described, Alday et al. designed a 1-dimensional quantum integrable system such 
that its free energy computes the value of the remainder function R^ 1 at strong coupling. In 
the 1-dimensional theory one can tune n — 5 complex mass parameters and the same number of 
real chemical potentials. The dependence of the free energy on these parameters captures the 
dependence of R^ n > on the relevant kinematic variables. 

Within the 1-dimensional quantum system it is natural to consider a limit in which the masses 
are sent to infinity or, equivalently, the volume R of the 1-dimensional space becomes large. In 
such a limit, all computations simplify once finite size corrections can be neglected. This applies 
in particular to the free energy of the ground state. The main goal of our work is to show that 
such a large volume limit of the 1-dimensional system possesses a nice re-interpretation in terms 
of the 4-dimensional gauge theory: It corresponds to the multi-Regge limit. Put differently, the 
map between 4-dimensional kinematic variables and parameters of the 1-dimensional system 
sends the multi-Regge regime to a point at which all the mass parameters become large. A 
more precise formulation of the limit in the 4-dimensional gauge theory will be given in section 
2. The identification (5.2) of the corresponding regime in the auxiliary quantum system is one of 
the main results of this work. It is derived in sections 4,5 and generalizes previous observations 
[24] for the case of six gluons to an arbitrary number of external particles. 
If we were only interested in the ground state energy of the system, the large mass limit would 
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be of limited interest. But it turns out that some excited states of the 1-dimensional quantum 
system also play an important role. In order to see them enter let us recall that the Regge limit 
of scattering amplitudes can be taken in different regions of the kinematic variables, such as the 
Euclidean region, the physical region where all energies are positive or 'mixed' physical regions 
with positive and negative energies. The limiting value of the remainder function depends on 
the region. In fact, when we pass from one region into another by continuation in the kinematic 
variables, the amplitude picks up Regge cut contributions that may have a non- vanishing high 
energy limit. In this sense, values of the remainder functions in the multi-Regge limit of differ- 
ent kinematic regions probe the analytical structure of the amplitude. One may wonder what 
all this corresponds to within the 1-dimensional auxiliary system. Since the kinematic variables 
are mapped to system parameters (masses and chemical potentials), we must vary the latter 
in order to move from one region of the kinematic variables to another. In the 1-dimensional 
system such a variation of system parameters can lead to a pair-wise creation of excitations 
above the ground state [25, 26]. The energy of such excited states may be non-zero in the large 
volume limit. Since the energy in the 1-dimensional system is related to the remainder function, 
excitations of the auxiliary model correspond to Regge cut contributions in the gauge theory. 
One example of this phenomenon was worked out in [24] for the case of six external gluons. 
Combining the insights from the previous two paragraphs we must address the challenge of 
computing excitation energies in the infinite volume limit. When finite volume corrections can 
be neglected, excitation energies are determined by a set of algebraic Bethe ansatz equations. 
These replace the more complicated non-linear integral equations that govern a 1-dimensional 
integrable system at finite volume. The data that enter the Bethe ansatz equations, namely the 
momenta and 2 1—7-2 scattering phases, can be derived from the non-linear integral equations. 
We will explain the general construction in section 6. In the case of six external gluons the 
derivation of the relevant Bethe ansatz is particularly simple so that we can make things very 
explicit. Starting from n = 7, an interesting new feature appears. In going to the multi-Regge 
regime of the 1-dimensional quantum system experiences wall-crossing, i.e. the associated non- 
linear integral equations pick up additional terms which we will compute in section 6. One can 
perform the large volume limit of such modified integral equations, but that leads to modifi- 
cations in the Bethe ansatz, as well. More explanations and explicit formulas are included in 
section 6 along with a sketch of how one may proceed to bring the Bethe ansatz equations for 
n>7 into the standard form. 

From the point of the auxiliary quantum integrable system, the multi-Regge limit is opposite 
to the high-temperature (small mass, R) limit of the Y-system that was considered by Alday 
et. al. [23] and then studied in more detail in [27-29]. In terms of the 4-dimensional kinematics, 
the high-temperature limit corresponds to the case where the gluon momenta pi form a regular 
polygon that can be embedded in a subspace R 1 ' of the full momentum space R 1 ' 3 . Another 
limiting regime of the kinematic variables is probed by the operator product expansions (OPE) 
of polygonal Wilson loops, see [30-32]. The information encoded in such Wilson loop OPEs 
seems more closely related to the multi-Regge limit, though the precise link is a bit difficult to 
establish even at weak coupling [33]. 
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2 Multi-Regge kinematics 



In this section we discuss the relevant variables and kinematics necessary for the description of 
2 — > n — 2 scattering in the multi-Regge limit. The multi-Regge limit is characterized by the 
behaviour of a particular set of Mandelstam invariants. The remainder function of scattering 
amplitudes in N = 4 SYM theory, on the other hand, depends only on very special cross ratios 
of Mandelstam variables which are invariant under dual conformal symmetry. Our task here is 
to describe the multi-Regge limit in terms of such cross ratios. 

2.1 Kinematic variables 

We are interested in the scattering of two incoming particles with momenta —p\ , —pi resulting 
in a (n — 2)-particle final state with outgoing momenta P3,.-.,p n as shown in figure 1. It will be 
convenient to label momenta pi by arbitrary integers i such that pi+ n = p%. 



2 




Pi Pn 



Figure 1. Kinematics of the scattering process 2 —> n — 2. On the right-hand side we show a graphical 
representation of the dual variables Xi. 

In the context of M = 4 SYM theory it is advantageous to pass to a set of dual variables X{ 
such that 

p i = x i - 1 -x i . (2.1) 

The variables Xi inherit their periodicity Xj+ n = Xi from the periodicity of the pi and momentum 
conservation. Let us also introduce the notation Xij = Xi — Xj. The Xij provide a large set of 
Lorentz invariants x\j = x^. 1 When expressed in terms of the momenta, these read 

x% = ( Pi+1 + ■ ■ ■ + Pj ) 2 . (2.2) 
1 Throughout this paper we use the metric (1, — 1, — 1, —1). 
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Obviously, only 3n — 10 Lorentz invariants are independent. Throughout this text we shall 
describe scattering processes through the 3n — 10 independent variables 



t r — x l r ^_2 i (2-3) 
s r = x r+1>r+3 , (2.4) 
r, p = ^-y'-^+i > (2 5) 

x p— 2,p+l 

where r = l,...,n — 3 extends over all i-channels and p = 4, ...,n — 1 labels the produced 
particles. In discussions of the multi-Regge limit it is actually quite common to use the cosines 
z r = cos 9 r of the scattering angle 9 r defined in the CM-system of the momentum q r instead of 
s r and to replace our variables r\ p by the so-called Toller angles uj pp+ i. With these choices, the 
multi-Regge limit is obtained by sending z r — > oo with t r and oo pp +i held fixed. The variables 
we defined in eqs. (2.4) and (2.5) are more convenient (for a pedagogical discussion see [34]). In 
the variables (2.3)-(2.5) the multi-Regge limit is taken by sending s r to infinity while keeping 
both t r and r/ p fixed. 

As we recalled in the introduction, the missing remainder functions for gluon scattering in N = 4 
SYM theory possess dual conformal symmetry. In other words, they only depend on conformally 
invariant combinations of Mandelstam variables. For an n-gluon scattering amplitude there are 
only 3n — 15 independent conformal invariants. We choose to work with the following set of 
cross ratios: 

2 2 

X (7+l,(7+5 X <7+2,<7+4 fn a , 

• i (X+2,<7+5 x o-+l,o-+4 

_ x <r+3,n x l, CT +2 (r> 7 , 

U2a — —9 o , 

x cr+2,n x l,o-+3 

X 2,(7+3 X 1,(7+4 , n oN 

X 2,(7+4 X l,(7+3 

where <r = 1, ...,n — 5. For the 2 — )• 5 scattering process, a convenient graphical representation 
of the cross ratios is displayed in figure 2. Our main task now is to analyze the behavior of 
these 3n — 15 cross ratios in the multi-Regge limit. 

2.2 Scattering in the center-of-mass system 

To study the behaviour of the cross ratios in the multi-Regge limit we need some preliminary 
results which we will obtain by specializing to the CM-frame, writing the results in Lorentz 
invariant form. In our analysis throughout this subsection we shall study the behaviour of the 
set of sub energies 

Si ...j = ( Pi + ■ ■ ■ + Pj ) 2 . (2.9) 

Note that these subenergies are the same as the Lorentz invariant variables xfj = Sj+i...j we 
introduced in the previous subsection. The only reason we change notation here is to give the 
equations in this section a more familiar form. Subenergies Sjj+i for two adjacent particles of 
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1 2 3 4 5 




Figure 2. Graphical representation for the cross ratios of the 7-point amplitude. 

momenta pi and pi + \ are related to our variables s r through s r = s r+ 2 r +3- Furthermore, the 
total energy s = (pi + P2) 2 of the process is given by s = s 3 ... n . 

In studying the Regge behaviour of the subenergies, it is useful to introduce the Sudakov 
parametrization 

Qr = S r pi + 7 r p2 + Qr±, (2.10) 

for r = 1, . . . , n — 3. Here, p\ and p2 are light-like reference vectors from which we define our 
incoming momenta as p\ = —pi, P2 = —p2- They obey 1p\p2 = ^P\P2 = s, and the transverse 
part, q r ±, is orthogonal to both p\ and p2, i.e. piq r ± = f>2Qr± = 0. A convenient frame is the 
CM-system of the incoming particles 1 and 2, with momenta p\ and p2 along the z-direction. 
We can determine the Sudakov parameters j r and 5 r by considering the following subenergies, 

S3-H-2 = (P3 H \~Pr+2) 2 = (~P2 ~ Ql + Ql - Q2 H h q r -l ~ q r ) 2 = (j>2 + qr) 2 

s r+3 ... n = (p r+3 h h Pn) 2 = {q r - q r +i + q r +i \- q n -3 - pi) 2 = (pi - qr) 2 ■ 

Using that qf = tj we find 

s 3 ... r+2 = 2q r p 2 + t r = -5 r s + t r =>■ 5 r = — — S3 "' r+2 , (2-11) 

s 

as well as 

s r+3 ... n = t r - 2q r pi => 7 r = . (2.12) 

s 

Up to now all the identities have been exact. Now we would like to continue considering the 
multi-Regge limit which, as defined in the previous subsection, amounts to sending all pairwise 
energies s\, . . . , s n -3 to infinity, while keeping both t- and rj- variables fixed. In this paper we 
will restrict ourselves to the physical kinematic region where all energies are positive and all t r 
negative (in a future study we we will consider also analytic continuations into other 'mixed' 
physical regions where some energies are negative). For the subenergies introduced in eq. (2.9) 
the multi-Regge limit implies 

S > S 3 ... re _i, S 4 ... n > S 3 ... n _ 2 , S5-TI > • • • > S >-ii,..., -t n _ 3 . (2.13) 
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From this well-known hierarchy of energy variables along with eqs. (2.11) and (2.12) one deduces 
the strong ordering of the Sudakov parameters 7$ and 5i 



1 > 71 > 72 > ••• > 7n-3 (2-14) 

and 

1 » -5 n - 3 » -£ n _ 4 > • • • > (2.15) 
As a simple consequence of eqs. (2.11) and (2.12) we note that in the multi-Regge limit 

t r = ql = s-jrSr + <2v 2 _l = q^ ± , (2.16) 

where we could drop the term s~f r 5 r = —s~ 1 S3... r+ 2S r +3-n because of the strong ordering (2.13) 
of subenergies in the multi-Regge limit. In conclusion, the finiteness of the t r in the multi- 
Regge limit implies that the transverse components of the q r stay finite. Since p r +3 = q r — Qr+i 
this is also true for the transverse components of the momenta p r +3 f° r r = 1, . . . , n — 4. We 
can compute this finite quantity from the mass-shell conditions of the produced particles with 
momenta p^, . . . ,p n -i- 

= P p = {q P -3 ~ Qp-2) = -sjp- 3 5 p -2 + P p ± = h P p ± 

^ Sjt ^ = -p 2 p± =Pp 2 . (2.17) 

We are now prepared to look at the subenergies in the multi-Regge limit. Let us begin with the 
subenergies formed by two adjacent particles. We have expressions of the form 

s r = (g r _i - q r +i) 2 = -S7 r _i5 r+ i + (p r+2 +p r+3 )i 



- + (p r+2 + Pr+3 )l (2.18) 



s 

where r runs over r = 2, ...,n — 4. Similarly, we can determine the leading terms in the multi- 
Regge limit of the subenergies for three adjacent particles, 

V-ipp+i = (q P -4 ~ Qp-i) 2 = -S7p-4<5 p -i + (Pp-i + Pp + Pp+i)± 

= S3 - p+lS s P - 1 - n +(Pp-i+Pp + P P+ i)l (2.19) 

for p = 5, n — 2. As an application of these results we can now express the variables r] p 
in the multi-Regge limit through the momenta of produced particles. In order to do so, we 
express the basic definition (2.5) of the r/ variables through the subenergies (2.9). All three 
subenergies that appear in the expression can then be replaced by their leading behavior in 
the multi-Regge limit, i.e. the first term in eqs. (2.18) and (2.19), respectively. Comparing the 
resulting expression with the result (2.17) we arrive at 

Vp = ^^ = -P 2 p± =p 2 p. (2.20) 

Sp— lpp+l 
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We can apply this result to derive a relation between our rj variables and the (azimuthal) Toller 
angles between between adjacent vectors q r , q r +\- The angle between vectors q r and q r +i can 
be determined by computing 



Pr+3 = (Qr ~ qr+lf = |*r| + |*r+l| ~ 2 sf \t r \\t r +i \ COS 9 r ,r+l (2.21) 

for r = 1, ... — 4. According to eq. (2.20), the quantity p r +3 coincides with the variable r] r +3 
in the multi-Regge limit. Hence, we obtain 

a \tr \ + l*r+l| _ Vr+3 /„ 00 n 

cosp rjr+ i = z . {2.22) 

2\/\t r \\tr+l\ 

Below we shall need an explicit expression for the sine of the (azimuthal) Toller angle Q r ,r+\ in 
the multi-Regge limit. It is given by 

• a A ~2a ~ Ml*r|) \tr+l\iVr+3) h 
sin 6^ r+1 = J 1 - cos^ 6 r , r+ ± = — , (2.23) 

2-\/|ir||*r+l| 

with 

A 2 (|£ r |, |t r +i|, rj r+ 3) = 2\t r \\t r+ i\ + 2r] r+ ^\t r \ + 27/ r +3|t r -|-i| — |t r | 2 — |t r +i| 2 — Vr+3- (2.24) 

From these angles between adjacent vectors q r and q r+ \ it is straightforward to compute the 
angles between arbitrary vectors q r and q r i due to the 2-dimensional kinematics in the multi- 
Regge limit. 

So far, we have only looked at subenergies for up to three particles. But it is clear how to 
continue the analysis. Generalizing the derivations of our eqs. (2.18) or (2.19) we obtain 

^ S r+ 2... n S3-r'+3 . I n2 / n qc;\ 

S r+ 2...r'+3 = V (Pr+2 + ■ ..p r '+3)± ■ {2.25) 

S 

Comparison of the leading terms allow us to conclude that 

S r+2 ... r '+3 = (Pr+2 H hp r ' +3 ) 2 = ~ ~ • (2-26) 

rir+3 ■ ■ ■ Vr>+2 

Note that for r' = r + 1, i.e. when the subenergy on the left hand side involves three particles, 
the relation is exact and not restricted to the multi-Regge limit. For more than three particles 
contributing to the subenergy, on the other hand, the result (2.26) only describes the leading 
term and it takes more effort to determine the subleading term from eq. (2.25). In the next 
subsection we only need to determine the subleading contribution for very special combinations 
of subenergies. We postpone further discussion of such subleading terms until we have spelled 
out the relevant combinations. 

2.3 Cross ratios in the multi-Regge limit 

Let us now look at the behavior of the cross ratios defined by eqs. (2.6)-(2.8) in the multi-Regge 
limit. Combining eq. (2.26) with eq. (2.2) we obtain the leading term in the multi-Regge limit 
of the basic Lorentz invariants xf +r 3+r , : 

x 1+ r3+r > = s r+2 ... r '+3 = • (2-27) 

Vr+3 ■ ■ ■ Vr'+2 
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If we insert this asymptotic behaviour into our definitions (2.7) and (2.8) for the cross ratios 
U2(j and u-sa we conclude 

u 2a = ~ 2 = 7 ' (2.2X) 

Icr+l X <j+2,n r cr+l S cr+1 
2 

t<r+2 X 2,a+3 ^ ta+2 Vo+3 , n nr .s 

u 3(7 = 2 — = ; • (2.2yj 

£<t+1 3-2 o+4 "<r+l ^0+1 

When we send s a to infinity to reach the multi-Regge regime, both sets of cross ratios go to 
zero. Since the coefficients of 1/ ' s a depend on the t— and rj variables only, the ratios 

U2a ^ (2 30) 

U3a tcr+2 "Ha+3 

approach a non- vanishing constant value in the multi-Regge limit. Functions of the cross ratios 
that remain finite in the multi-Regge limit can therefore depend on the ratios U2a/u 3{T . 
Let us now look at the remaining set of cross ratios u\ a . Once more, we can rewrite our definition 
(2.6) in terms of the energy variables (2.9) using eq. (2.2) to obtain 

u\ a = . (2.31) 

■V+3- ■ -cr+5 s <r+2- ■ -0-+4 

If we now insert the limiting behavior (2.26) for all four subenergies we see that all variables 
u\ a behave as u\ a = 1 + . . . . In order to obtain the leading non-trivial term in the multi-Regge 
limit we must work a little harder. 

To this end we write, in analogy with eqs. (2.18)-(2.19), expressions for the subenergies s CT +3 o+4 5 
s CT + 3 ... CT+5 , s CT+ 2... CT+ 4, and s a+2 ... a +5- Beginning with the relation for s a+3 CT+4 

1 = *:<■■■» ,;{...„ ^ (Pa+3 + Pa+4:)± ^ ^ 

SS a +3 o+4 s o+3 cr+4 

we write 

_ S3...cr+4S(T+2-n S3---o+5So+3---n SS a+ 2--- a +5 S a+ 2-g+4Scr+3-cT+5 
SS a +2-(7+4: SS CT +3-o+5 S 3 -a+5Sa+2-n «o+3 <7+4S(7+2-<7+5 

+ Cp ff +3+P ff +4)l j {2 33) 

s cr+3 cr+4 

On the right-hand side, the fourth fraction equals , and the first three fractions are equal 
to unity, arising from the equations for s CT +2... a+4, s<r+3— <r+5> and s a+ 2-a+5, with corrections of 
the order C(l/s - + 2... -+4), 0{lj s^z—a+b), C(l/s - + 2---o-+5), respectively. Compared to the last 
term in eq. (2.33), these corrections can be neglected and we are left with 

s a+3 cr+4 

with further corrections being smaller than 0{1/ s a+3a+ ^). The numerator of the correction 
term in eq. (2.34) can be expressed in terms of Lorentz invariants. To do so, we use the results 
of section 2.2 by writing 



(p a+3 + p a+i ) 2 = (q a - q a +2) 2 = \ta\ + IW2I - 2^1^11^+21 cos(6» CTj(T+ i + 6» CT+ i iCT+2 ), (2.35) 
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which finally gives (p a +3 + P0-+4) 2 = Pa with a set of functions of the t and rj variables defined 
by 

pa{t,V) ■= \U\ + \ta+2\ ~ 2 y/ \t a \ \t a+2 | = = 

\ 2^/|r cr ||r (T _)_i| 2y|t - + i||t cr _|_2| 



A(|to- 1 , 7/0-4.3) A(|tg- + i|, 1 ? 7 ?o-+4) 

2v / |*cr||i ( T+i| 2^/|to- + i I |t cr _ l _2 1 



(2.36) 



We can now summarize the findings of our analysis on the multi-Regge limit of the cross ratios 
u\ a through 

u la -l = p a (t, T])/Sa+1 (2.37) 

where we also changed notations back using s a+ \ = 50-4.30-4.4, as mentioned before. As in the 
case of the cross ratios u 2a and the leading correction to u\ a — 1 vanishes in the multi-Regge 
limit. But the following ratios remain finite 

gig - 1 ^ pJhvytg+i ^ 2 38 ^ 

U2a taVa+A 

This concludes our description of the kinematics in the multi-Regge limit. 



(2.39) 



3 The n-gluon thermodynamic bubble ansatz 

The main goal of this section is to review the Y-system for the computation of n-gluon ampli- 
tudes at strong coupling [22, 30]. In the first subsection we explain how the most interesting 
contribution to the scattering amplitude can be computed by solving a system of non-linear 
integral equations (NLIE). Then we relate the parameters of the NLIE to the cross ratios that 
were introduced in eqs. (2.6)-(2.8). 

3.1 Amplitudes and the Y-system 

We are interested in the calculation of scattering amplitudes in N = 4 SYM at strong coupling. 
To leading order they are given by 

Amplitude ~ e~ 2^ A , (3.1) 

where A is the area of a minimal surface in AdS§ with piece-wise light-like boundary. A general 
prescription for the calculation of this area A is given in [23, 35]. It contains a number of 
different pieces, including a divergent BDS-like term and a number of finite contributions. All 
but one of these terms can be spelled out explicitly. The remaining one is also known, but 
it is characterized somewhat indirectly through the solution of a coupled system of non-linear 
integral equations. Because of the resemblance with the way one describes the free energy of 
a 2-dimensional quantum integrable system, this contribution to the area A has been dubbed 
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^4f ree . In our analysis of the multi-Regge limit we can restrict to the discussion of this free 
energy contribution since the remaining terms are straightforward to include. 
For our study of scattering amplitudes in the multi-Regge regime we need some more background 
on j4f ree - It can be calculated from a set of functions Y ajS with a = 1,2,3 and s = 1, . . . , n — 5, 
which are determined as solutions of the following set of integral equations: 

logYi jS = -m s cosh - C s - -K 2 * /3 S - K\* a s - ^K 3 * j S) (3.2) 

logY2 j s = — m s V2cosh0 — K 2 * a s — K\ * (3 S , (3-3) 

logY 3jS = -m s cosh 9 + C s - -K 2 * /3 S - K x * a s + ^K 3 * 7 S , (3.4) 

where K * f denotes the convolution integral 

CO 

J de'K(o-e')f(o') (3.5) 

— oo 

and a s , f3 s , ^ s are given by 

, (1+Y 1|a )(l+Y 3 , a ) 
a * = lQ g 7i I v »i l v \) (3-6) 



'(1 + Y 2>a _i)(l+Y 2) , +1 )' 



\2 



P> 0g (l + Y liS _ 1 )(l+Y M+1 )(l + Y3, s _i)(l+Y 3 , s+1 )' 1 • ' 

^ (l + Y 1 , a _ 1 )(l+Y 3|8+ i) 
7s = l0g (l + Y 1 , s+1 )(l + Y 3 ,_ 1 )- (3 - 8) 

The kernel function lf a are known to take the form 

K - 1 1 if - ^ C ° Sh 9 K - 1 t h 29 (3 9) 

2-7T COsh 9 ' 7T COsh 29 ' 7T 

Furthermore, m s and C s are constants that we need to determine in the following. These 
equations can be used for |Im0| < j. For larger values of 9, we can either pick up pole 
contributions from the kernels or use the recursion relation 

i+Y^nfi+Y^ 



■ a,s+l M x ' x 4-o.s— 1 , 
= [r V ±2] / \ K x 7 

Yf-a 2 ! ( 1 + ^F±W ) f 1 + \7fW 



a + l,s / \ a — 1, 



where we introduced the symbol Ya} s (9) = Y ajS (9 + irir/A) for Y- functions with arguments 
shifted by multiples of wr/4. For the moment let us consider the m s as complex parameters 
while we take C s to be real. Consequently, the total number of real parameters in the eqs. 
(3.2)-(3.4) is 3(n — 5), matching the number of independent cross ratios for n-gluon scattering. 
The precise relation between the cross ratios and the parameters m s , C s will be addressed in 
the next subsection. Right now it suffices to keep in mind that the parameters m s = m s (u aa ) 
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and C s = C s (u aa ) in the Y-system need to be adjusted as we vary the kinematic variables. 
Let us add a few more comments on complex mass parameters m s . The correct way to interpret 
the Y-system in the presence of complex masses is through the following substitutions in the 
original equations: 

m s ^\m s \, Y a , s (0)^Y a!S (0 + ^ s ), K?${0 - 9') (9 - 9' + i(<j> s - fo)). (3.11) 

Here we have split each complex m s into the real parameters \m s \ and the phase (f) s . For 
1 0s — 0s' I = n • \ the kernels become singular, and we have to pick up the corresponding poles. 
Such large phases are going to play an important role later on. 

Once the solution Y a ^ s of the Y system has been found, we can compute the quantity ^4f ree 
through the simple prescription 

A ixcc = Y, \ ^|m s |coshMog[(l+Y M )(l + Y 3 , s )(l + Y 2iS ) v/1 ] (9 + %$.). (3.12) 

Note that A{ ree depends on the kinematic variables u aa of the scattering process through the 
parameters m s , 4> s and C s of the auxiliary quantum integrable system. 

3.2 Y-system and cross ratios 

In this section we relate the cross ratios to the value of the Y-functions at special values of the 
spectral parameter 9. To do so, we follow [23] and define 

1 



U[r] := i + l = i + 



1 + 



(3.13) 



=mit/ 4 Y 2 ,s 

for s = 1, ... ,n — 5 and r any integer. These quantities possess a rather simple relation with 
the cross ratios. When both indices of U are even, one has 

TT [2p] _ x -k+p,k+p x -k+p-l,k+p-l fo-\A\ 
U 2k-2 - 2 -2 ■ \ 6 - lli ) 

-fc+p-l,fc+p -fc+p,fc+p-l 

Here 2k — 2 is the number of cusps between the pairs {x-k+ p -i, x-k+ p } and {xi~+ p -i, xj-+ p }, 
see [23] for details. For sites separated by an odd number 2k — 1 of cusps, the relevant relation 
reads 2 

TT [2p+l] _ X -k+p,k+p+l X -k+p-l,k+p /„ , ^ 

U 2k-1 ~ Z2 12 ' ^- iD J 

-fc+p-l,fc+p+l -fc+p,fc+P 

We can now insert the general relations (3.14) and (3.15) into our definition of the special cross 
ratios eqs. (2.6)-(2.8) to obtain 

_ X a+l,(T+b X <T+2,a+A _ I ,+2(7+7]^ _1 _ 1 2,1 
■ i '(7+2, ( T+5 X o-+l,o-+4 v / 1 + r 2 L 1 

2 2 y-[°'+ 4 ] 

X cr+2,n X l,o-+3 X 7 1 + Y 2 L j(T 

V+6l\ _ X 2,a 



U2o 



= I = 7TL< 

T 2 T 2 \ J 1 , 

X l,o-+3 X 2,o-+4 v 7 1 + Y 2 L 
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The Y- functions appearing in our theory depend on the 3n — 15 parameters m s ,C s ,(j) s in the 
non-linear integral equations. Hence, Y a ' s which are just shifted Y-function evaluated at the 
origin of the 9 plane, are functions of these parameters. The equations (3.16) describe the 
transformation between the parameters in the non-linear integral equations and the cross ratios 
of the scattering process. We can invert them, at least numerically, to determine m s , C s and 4> s 
from the kinematic invariants of the n-gluon system. 

The formulae derived above for the cross ratios give rise to large values in the upper index of 
the Y-functions. As we shall see later, this is a bit of a nuisance for practical computations. 
It is therefore useful to observe that there exist two symmetries which may be used to reduce 
the value of the upper index. These symmetries have their origin in the Z4-symmetry of the 
underlying Hitchin system. The first of these symmetries reads 

Uf = Ul l±2n] . (3.17) 

Note that such a symmetry must necessarily hold in order for the identification (3.14) and (3.15) 
with cross ratios to be consistent with the symmetry xi +n = x% of the x-variables. A second 
useful symmetry of the quantities is given by 

D? = Ultts ■ (3-18) 

In this case, we must accompany the shift in the upper by a reflection in the lower index. Once 
again, the corresponding symmetry of cross ratios is easy to verify. With these symmetries, it 
is always possible to reduce the absolute value of the upper index of the Y-functions to I S I or 
lower. In order to achieve further reduction, one can employ the recursion relations (3.10). 

4 Multi-Regge limit of the TBA for n < 7 gluons 

The multi-Regge limit was defined in section 2 through the dynamical invariants of the scat- 
tering process as a limit in which the s-variables are sent to infinity while t- and r/-variables 
are held fixed. We have also analyzed how the special cross ratios (2.6)-(2.8) behave in the 
limit. In section 3 we then went on to discuss the relation (3.16) between cross ratios and the 
parameters of the non-linear integral equations. Our next task is to understand which limit of 
the parameters m s , 4> s and C s has to be taken in order for the cross ratios to show multi-Regge 
behavior. The case with n = 6 has been treated before [24] and is relatively simple to analyze. 
We will review some formulas in the next subsection before turning to n > 6 gluons. Beyond 
n = 6 there are some important new features in taking the multi-Regge limit. We will explain 
these first for the example of n = 7 before we delve into a general analysis in the subsequent 
section. 

4.1 Review of the hexagon n = 6 

In order to understand the basic steps of our analysis we would like to review briefly how things 
work in the case of 6 gluons [24]. For six points, there are three independent cross ratios: 

Jl Jl ^[-3] 2 Jl y-[-l] Jl Jl yM 

u ^ _ x 35 x 26 _ x 2 U2 _ x 13 x 46 _ x 2 ^ _ x 15 x 24 _ J 2 

x lA i + y 2 M] ' 2 x \A i + y 2 M] ' 3 x \A i + yJ 11 ' 
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We have omitted all the second indices for both u and Y because a = 1 is the only value it can 
take when n — 5 = 1. From the kinematic analysis in section 2 we know that U\ — > 1 in the 
multi-Regge limit, while U2 and tend to zero. Comparing the expressions for the Y-functions 
with this limit, the authors of [24] show that 

m — > oo, <]) — > 0, C = const. (4-1) 

is the appropriate limit one has to perform in the non-linear integral equations in order for the 
cross ratios to assume their limiting values in the Regge regime. When the limit (4.1) is taken, 
the integrals in the NLIE (3.2)-(3.4) may be neglected. Consequently, we obtain a set of explicit 
expressions for the form of the Y-functions in the limiting regime: 

dd g — m cosh(0— «</>)— C Y2($) ~ g— mv / 2cosh(0— «0) ~Y^(0^ — e~ m cos ' 1 (^ — *</ > )+C 

Recall that these expressions should only be used for |Im0| < tt/4. Outside this fundamental 
strip one needs to apply the recursion relations (3.10) to bring the arguments back into the 
strip. Before we insert these expressions into the formulas (3.16) let us replace m and <j) by the 
new variables 

e = e -mra *, w = e msin(t ', (4.2) 

which behave as e — > and w —> const, in the limit (4.1). In terms of these new parameters, 
the cross ratios (3.16) can be expanded as 

Ul = 1 - ( w+ — + 2 cosh C ] e + C(e 2 ), u 2 = we + C(e 2 ), u 3 = - + 0{e 2 ). 

\ w J w 

While U2 and only involve the Y function Y2 in the fundamental strip, we need to use the 
recursion relation (3.10) to find u±. Hence, all the three cross ratios indeed show their Regge 
behavior (2.28), (2.29) and (2.37). 

4.2 Multi-Regge limit for n = 7 gluons 

For the 7-point amplitude, the cross ratios are given in terms of the Y-functions as 

jr[2] y[-2] y[0] 

2,2 2,2 2,2 ( A Q \ 

Ull = 7~Wv U2l = T7W^' 31 = iT7m' (43) 

vi-3] v m 

2- 1 2* 1 2' 1 (A A\ 

U\2 = p^T, U 2 2 = TTTi «32 = [TT • (4-4) 

i + Yl 1 3] i + Yl, 1] i + yPj 

We are going to demonstrate that that the cross ratios obtained from the Y-system for n = 7 
display multi-Regge behavior if we take 

m s — > oo , C s = const. 

01^0 , ^2^-\- (4.5) 

Note that the limiting value for the second angle fa is non- vanishing. In the subsequent section 
we shall argue that in this limit, all integral contributions may be neglected and that the 
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prescription (4.5) is the only one that provides the correct Regge asymptotics of cross ratios. 
For the moment let us just check how things work with the limit we propose. As in the discussion 
for n = 6 we shall switch from variables m s and 4>s to 

Cl = e~ mi cos ^ , Wl = e mi sin ^ , (4.6) 

e 2 = e ~ m2 cos (f +^) , W2 = e m2 sin (f +^) . (4.7) 

In the above mentioned limit (4.5), these quantities behave as ej — > 0, Wi — > const. Let us begin 
our study of cross rations the simplest cases, namely the two cross ratios 

^ " ITWJ = 1 + e"^ -(i"*0 " ^ " ^ ^ } ' ( } 

u 2 2 = Ar^r = — 7= 7-, 7 = ; = ei^i + 0(e 2 )- (4.9) 

Up to this point, things work pretty much the same way as for the cross ratios 112 and U3 in the 
case of n = 6. The next cross ratio we want to look at is 

y[°l -v^macos^ £L 

U31 = ?i2 ^ _J = = 11 + Ote 2 ) (4.10) 

which is the first one to contain <j)2- It is this last computation that suggests for the first time to 
set the limiting value of 4>2 to fa = 7r/4. If we had set 02 = 0, for example, we would have been 
forced to omit the shifts by 7r/4 in the arguments of the trigonometric functions in eq. (4.7) in 
order to ensure finiteness of W2- Without the shifts, the Regge limit of U31 would have been 
given by €2- But since the definition of W2 and 62 included a shift, we had to add and subtract 
the limiting value 7r/4 of the angle 4>2 in the argument of the cosine. This is how we obtained 
the familiar looking Regge asymptotics of U31 even though the construction of the cross ratio 
only involved Y 2 °l with a vanishing upper index. 

All remaining cross ratios involve values of Y-functions outside the fundamental strip so that 
we need to make repeated use of the recursion relation. We see that 

1 4- 1 , 

and therefore 



2,2 \ * ^ yl-l] \ ^ 

J 3,2 J \ 1 1,2 



Y 



-2] 



-2,2 ^ £2^2 . rn( 2\ (Al-i\ 

U2i = T^T - l~, =e 2 w 2 + 0{e z ). (4.11) 

1 + Y 2 [ ~ 2 l 1 + e 2W 2 



Analogously, we obtain 



2,2 



^{l + W2 e~C2)[i + W2e c 2 y 



(4.12) 



- 16 - 



from which we find that 



Y [2] ( l 

«ll = ^^1- lw 2 + — + 2coshC 2 ) e 2 + 0(e 2 ). (4.13) 



i + y 2 [2 J v w-2 

[—31 

The last cross ratio we need is given in terms of Yi, t . The Y-functions appearing in the 
recursion relation are given below (without the full calculation). 

v [-2] ^ l + e c H 2 [-2] ^ l + e- c *e 2 



This then leads to 

z 2,l 



3] _ 1 + Y 2,2^ ^ (1 + e 2 w 2 ) 




and finally 



F [-3] , 

«12 = ^r^T - 1 - ( ™i + — + 2coshCi ) ei + 0(e 2 ). (4.14) 

i + r^ 31 V ™i 



This shows that the choice of parameters indeed gives the right limits for the cross ratios. 
4.3 Finding the correct limit 

Our analysis in the previous section was based on the claim that integral contributions to the 
NLIE can be neglected in the limit (4.5). Under this assumption we showed that our cross 
ratios display multi-Regge behavior. On the other hand we also suggested that the limit (4.5) 
was uniquely fixed by these requirements. Let us now discuss these claims in a bit more detail. 
It is clear from the relations (4.3) and (4.4) that for a cross ratio vanishing in the multi-Regge 
limit, the corresponding value of the Y-function has to vanish, as well. Values of Y-functions 
that appear in the cross ratios u\ a , on the other hand, must diverge in the multi-Regge limit. 
For example, from the limiting behavior U31 — > we conclude y\ ?\ — > 0. Therefore, log has 

to be large and negative. The same follows for the value logY^i^ considering that U32 — > in 
the multi-Regge limit. Because of the way the parameter m s enters into eq. (3.3) we enforce the 
desired behavior if we make m s very large. In this limit, the integrals appearing in eqs. (3.2)- 
(3.4) can be neglected because, assuming real masses for the moment, their leading contribution 
can be written schematically as 

J d9'K(9-9')log(l + Y(0')) J d6'K(0-9')log(l + e- msCOshe '^ -»■ 0, (4.15) 

as m s — > 00. This remains true even after introducing complex masses, as is shown in section 
5.1. In this limit, we can now study the behaviour of the Y-functions and find constraints on 
the d>„ and C=. 
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For further illustration we now wish to look at the cross ratios u\\ and analyze how we ensure 
u\\ —> 1 in the multi-Regge limit. The relevant value of the Y- function is given by 

1 + e -V2micos(f-<£i) 

Y 2 2 - 7 —i \ w 7 s T -> °°. ( 4 - 16 ) 

e -v/2m 2 cos0 2 ( 1 -)_ e m 2 cos (f _ ^2)-C 2 \ fy _|_ e m2Cos(f-</> 2 )+C2 



If we want the Y-function to diverge, every term of the denominator has to vanish in the 
multi-Regge limit. Writing the denominator as 

e -V2m 2 co B <j> 2 + e v / 2m 2 sin^ 2 + 2 CO sh C 2 e~ m ' 2 COS (^ +02 ) (4.17) 

we see that 

<h€(-^)n(-n,0)n[- T ^ = (~,0), (4.18) 

at least if we assume that it will stay in the interval between — ir < fa < it. In order to find 
the specific value that fa should assume, we recall from eq. (2.30) that the ratio 1121/1*31 must 
remain constant in the multi-Regge limit. According to our equations (4.3), this requires 

» const. (4.19) 



y[°] e -2v / 2m 2 cos<^ 2 M _|_ e ™2 cos(f +<fe)-C 2 \ f-y j_ g m 2 cos(f +</> 2 )+C 2 



2,2 

Here we have also used that both F 2 2 2 ' and Kj°2 vanish in the multi-Regge limit. For the ratio 
of these values to be constant, at least one term of the denominator has to go to a constant 
with the other terms in the denominator going to zero. This suggests the two possible limits 
fa — > —j or fa —> \. The latter limiting behavior, however, is excluded by the constraint 
we derived from u\\ — > 1. A similar analysis is carried out for the remaining Y- functions in 
appendix A. There we derive all constraints one can put on the parameters of the NLIE. These 
are solved by fa — > and fa — > — j, as we anticipated above. And indeed in section 4 we saw 
that such a limiting behavior of the phases <j) s produces the correct multi-Regge behavior for 
all the cross ratios. 

The most surprising outcome of our discussion for n = 7 is that the limiting values of the 
phases need no longer be zero, in contrast to what we found for n = 6. One may generalize the 
arguments outlined here to larger number of gluons n > 8 to find that cp s = — (s — l)j appear to 
be the correct values for the phases in the multi-Regge limit. A full proof of this fact, however, 
requires a closer look at the Y-system for non-zero phases and possible residues. 



5 Multi-Regge limit of the TBA - general case 

As we have argued in the previous section, the phases fa may approach large values in the multi- 
Regge regime. Since the original Y-system is only valid as long as phases satisfy \ fa — fai\ < 7r/4 
we must be prepared to include additional contributions that arise from poles in the kernel, see 
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our comments in section 3.1 and [30]. In the following analysis we will make the ansatz 



m s — > oo, (5-1) 

<Ps = (5-2) 

C s = const. (5-3) 

for the multi-Regge limit of the parameters m s , <j) s and C s . Our main goal is to show that in this 
limit, the cross ratios (2.6)-(2.8) have multi-Regge behavior, i.e. that they behave as described 
in eqs. (2.28), (2.29) and (2.37). The contributions from poles of the kernel functions turn out 
to play a vital role in this analysis. Therefore, we begin with a few general comments on the 
behavior of the Y-system for large masses in the presence of large phases <p s . These then allow 
us to find explicit values for the physical cross ratios for an arbitrary number of external gluons. 

5.1 Large phase residues and multi-Regge limit 

In this section we revisit the structure of the Y-system equations in the presence of large phases 
4> s . Recall that in the presence of complex masses the equations can be written as 

logY 2jS (0) = -V2\m s \ cosh# + [ de'K 2 s 'f,(e -6' + i<f> 8 - ^)log(l + Y oV (0')), (5.4) 

a',s' •* 

where Y QjS (#) = Y a ^ s (6+i(p s ), cf. [23]. We only displayed the equations for Y2 jS here, but similar 
equations obviously hold for the other Y-functions, as well. Eqs. (5.4) determine the functions 
Y in the fundamental strip |Im#| < 7r/4 as long as the phase differences satisfy \cp s — <f> s i\ < 7r/4 
. For Im(^ + i<f) a — i(f> s ') = k ■ j, the kernels are singular and we have to pick up residues from 
these poles 2 . More precisely, the kernels K2, K3 become singular for k = In + 1, while K\ 
is singular for k = 2(2n + 1). In other words, at least one of the kernel functions has a pole 
whenever 

k e X := {2Z + 1} U {2(2Z + 1)} . 
Once we include the residues of the pole contributions, the integral equations are given by 

logY 2)S (0) = -\/2|m s | cosh 6 + ^ n^log (l -f Y avjSv (o + i(j> s - i<f) Sv - ^7)) 

V 

+ [ de ' K2 si (0-0' + i(j> s - # s ')log(l + Y ttV (6 1 ) ) (5.5) 

a>,s> J 

In eq. (5.5), it is understood that the sum in the first line extends over all the relevant pole 
contributions. We will discuss this in more detail below. As they stand, equations (5.5) are 
valid for arbitrary 9 in the complex plane. Of course, the number of terms in the summation 
over v depends on the value of lm(0). 

Before we proceed to analyze the behavior of eqs. (5.5) in the limit of large masses, let us have 



2 Note that we have to pick up a pole only if the kernel singularity is crossed. For values on the singularities, 
an iO-prescription can be used (cf. [23]). 
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a closer look at the Y-system for real argument 9. Recall from section 3.1 that the kernels K^'*, 
are non-zero only for s f = s or s f = s ± 1. Hence, the differences <p s — <j) s > that appear in our 
limit (5.2) are restricted to values \<f> a — (p s /\ < 7r/4. We conclude that, as long as 9 is real, no 
poles of the kernel functions are actually crossed in the limit (5.2) and hence the Y-functions 
Y Q S (0),0 G ffi, obey the original Y-system without additional residue terms. 
This comment becomes important in passing to the large mass limit of eqs. (5.5). Note that the 
integral in the second line extends over 9' along the real line. Hence, to evaluate the integral, 
we only need to know Y n / s i(9') for real 9' where it is unaltered by pole contributions and 
consequently still given by 

Y 2)S /(0') g-v^mslcoshe' ; Y a , ;S ,(0') e -\™ a \co S h8'+(a'-2)C s foro ' = lj3i 

in the limit of large masses \m s \ and for 9' £ IR. It follows once again that all integral contribu- 
tions to eqs. (5.5) can be neglected in the limit (5.2). Shifting back to the Y-functions, we see 
that after neglecting the integrals the equations look like 3 

logY 2 , s (0) = -v / 2|?n s |cosh(^-i ( / )s ) + ^n,log(l+Y a ^ s , y (tf-^-)) . (5.6) 

It is important now to describe the sum over pole contributions in some more detail. To begin 
with all labels only run over nine possible values. While a u is free to assume any of 

its three values, s v cannot deviate from s by more than one unit, i.e. 

a u = 1,2, 3 , s u E {s,s ± 1} . 

Given any such choice of {a v ,s u ) we compute the quantity k Si/ = ^Im(^ — i<ps v )- It determines 
the possible values of the integer 

Given (a u ,s u ) and k u we must still find the integer n v . It is determined by the form of the 
functions a s ,(3 s and 7 S in eqs. (3.6)-(3.8), the coefficients of these functions in the Y-system 
(3.2)-(3.4) and sign of k u . We will discuss some examples in the next subsection. 

5.2 Multi-Regge limit for n = 8 gluons 

As we shall discuss at the end of the next subsection, pole contributions from the kernels only 
start to enter the limit (5.2) starting from n = 8 external gluons. It is useful to analyze the 
n = 8 case first to illustrate the general formula (5.6). Generalizing the variables we used in 
our analysis of the multi-Regge limit for 7 gluons we introduce 

e s = e -m.oo B ((«-i)f+^) j ( 5 _ 7 ) 

Ws =e m.rin((.-l)f+^) > (5.8) 

where s = 1, 2, 3. In the limit (5.2) these variables behave as e s — > and w s —> const. 

[—21 

As an example, let us analyze the limit of F 2 2 . Since we want to determine Y 2i 2 at 9 = —in/ 2 
3 Note that the variable 9 was shifted by — i<j) s when we passed from Y a ,s to Y a ,s- 
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Figure 3. Residue structure for Y-functions in 8-point case. 



and <t> Sv = —{s u — l)7r/4, the parameter k Sv defined at the end of the previous subsection reads 
K s„ = s u — 3. In order for the interval (0, s u — 3) to have a non-vanishing intersection with X , 
we must have s v = 1. In the case at hand, the intersection consists of a single point k = — 1. 
Looking back at eq. (3.3) we note that only K2 possesses a pole at —iir/A. The coefficient of K% 
in the equation for log Y2,2 is —02- The latter contains only one contribution with 1 = s v = s— 1, 
namely the term — log(l + Y2 t i). Consequently, we find 



Note that in this particular case the correction term becomes trivial in the limit (5.2) since 
3^2 1^ = exp(— m s ). Analogously, we can analyze the residue structure of the remaining Y- 
functions. The result is shown in figure 3 which we explain in the following. Every node in the 
diagram represents the 3 values Y^ with a = 1,2,3 while s and k are kept fixed. Encircled 
nodes correspond to Y-functions that receive no corrections from residues. Arrows points to the 
nodes from which a given Y-function receives residue terms. Our result (5.9), for example, is 
represented by the arrow that connects the square box around Y^ 2 2 ' with the circle around Y^ ^ . 
Values of Y a , s that are not included in the figure not only receive residues with k v = ±1, but 
include higher values of k v . Those have a more complicated structure and are better bypassed 
through the use of recursion relations, if possible. We see that with the help of the Y-functions 
at 9 = —in/A we can determine all Y-functions at 8 = 0. Once we have constructed Ya jS ^ and 
Y^a} s , all remaining values can be reconstructed with the help of the recursion relation (3.10). 
A complete analysis gives the following results for the cross ratios in the limit (5.2), 

uxa = 1 - \ Wn-A-o H V 2 cosh C„_ 4 _ CT ] e„_ 4 _ (T , (5.10) 

U2a = ^n-A-aWn-i^a, (5-H) 
U 3a = e n ~4~aW~\^ a . (5.12) 
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Figure 4. Residue structure for n points. Y-functions not shown in the hgure correspond to gray boxes. 



So, we find that all cross ratios have multi-Regge asymptotics (2.28), (2.29) and (2.37). 
5.3 The general case of n-gluon scattering 

After discussing 7- and 8-gluon amplitudes we finally proceed to the general case. Our goal is 
to show that the cross ratios (2.6)-(2.8) that are obtained from the Y-functions through eqs. 
(3.16) show multi-Regge behaviour (2.28), (2.29) and (2.37) in the limit (5.2). As before, we 
can study the residue structure of the Y-system for phases 4> s = —{s — l)7r/4. The results are 
encoded in figure 4. Gray nodes correspond to Y-functions with a more complicated residue 
structure. We explain in appendix B how our pattern can be used to read off the resulting 
contributions. If a cross ratio (3.16) is obtained from the Y-function Y^j,, we have indicated 
this by putting the cross ratio into the node instead of the Y-function. We have already used 
the shift symmetry to put some of the u\ a in the first row. 

We see that the cross ratios U2 a and U3a lie ° n a diagonal whose 'residue flow^ stays within 
the diagonal. The endpoints of the residue flow are given by U31 and U2 n -5, which are free of 
residues and can be calculated directly to to give 

U31 = e n -5Wn-5i ( 5 - 13 ) 
U2n~5 = ei'^i- (5.14) 

It is not difficult to determine the remaining cross ratios U2 a and u^ a . Recall from eq. (3.16) 
that the cross ratios U2 a are given by 

y[-(n-A-a)] 

U ^ = ] ~ y [-(n-4-„)] • ( 5 - 15 ) 
1 I 2,n-i-a 
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According to our general expression (5.6), the values of the Y-functions that enter these cross 
ratios are given by 

Yl~ s] = e -v / 2m s cos(f+( s -i)f+^) . ( Residue terms) = e s w s ■ (Residue terms) . (5.16) 

We can conclude by iteration that the residue terms along this diagonal only introduce correc- 
tions of the order e 2 or higher and can be neglected. Hence, we obtain 

U2a = e n -4-aW n -4- a . (5-17) 

One can go through the same arguments to find 

u 3a = e n _ 4 _ CT w;^ 4 _ CT . (5.18) 

This finally leaves us with the cross ratios u\ a for which we have to work a little bit harder, as 
these variables all have a more complicated residue structure. We will make use of an inductive 
argument. Going from a n-point amplitude to a (n + l)-point amplitude introduces another 
row in our above pattern. However, the old Y-functions do not change their values as long as 
their residues are not affected by the new row. What changes, however, is the location of the 
cross ratios in the pattern. The most important change for our purposes is that all u\ a in the 
first row are shifted by two boxes to the left in our pattern, and a new cross ratio n 1 ^ n+1 )_ 5 ) 

\— 3l 

appears which is related to Y2 ± . This means, for example, that Ui^ n+ i^_^ in the (n + l)-point 
amplitude will take the value of «i n -5 in the n-point amplitude, since all residue are unaffected 
by the new row. Following our discussion in appendix B, it is clear that a Y-function will 
keep its n-point value if the diagonal from its node to the lower right hits the (n + l)-point 
U2o--diagonal at the position of U22, or above. It turns out that the last element, i.e. the element 
with smallest value of 9, in the first row that keeps its value is given by U13. This means, that 
all u\ a for a = 3, ...,n — 4 can be calculated using the known n-point values. However, u\\ 
and u\2 are then given by the symmetry (3.18) relating u\ a U\i n -^-a\- We have established 
above that 

ui a = 1 - I w n _ 4 _ CT H h 2coshC n _ 4 _ (T ] e n _ 4 _ CT (5.19) 

for the 7- and 8-point amplitude. By the previous argument, the same remains true for all 
values of n. This shows that the n-point solution is a (n — 5)-fold copy of the 6 point solution. 
In conclusion, we now confirmed our identification of the multi-Regge regime with the limit 
(5.2) of the parameters in the Y-system. 

5.4 The 7-point case revisited 

In the light of the general result derived in the last section, let us revisit the 7-point case to study 
explicitly the equations governing the Y-system in the fundamental strip and to understand 
why no large residues had to be considered in our earlier analysis. Recall that for 7-points the 
multi-Regge limit corresponds to taking m s large, C s constant, <f>i = and fa = — vr/4. For 
— 7r/4 < Im# < we have that Im(# — icp s ) G (— 7r/4, vr/4). Since none of these terms crosses a 
multiple of 7r/4, no singularities arise and we do not have to pick up any residues. Therefore, the 
Y-system (5.4) is valid without any modifications. For < Im# < 7r/4, however, Im (9 — ifa) 
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crosses 7r/4 and we have to pick up residues. Following the procedure outlined in the previous 
subsection, we find the following equations: 

logYn (0) = - |mi| cosh (0 - ifa) + Y, I d0 ' K i,a' i d ~ d> ~ lo S i 1 + Y *V' i 9 ' + *&0) 

a',s' ^ 

-C 1 + log(l + Y 12 (e-<^)), 
logY 21 (0) = - v^|mi| cosh (0 - »^) + J] / dfl'ifjj; (0 - 0' - log (1 + Y a ,, s , (0' + i0 fl ,)) 

a', a' 

+ log(l + Y 22 (0-ij)), 
logYsi (0) = - \mi\ cosh (0 - ifa) + / ^'^S' " " lo S ( x + Y a',*' + i^a')) 

a',s' J 

+ C 1 + log(l + Y 32 

logY 12 (0) = - |m 2 | cosh (0 - i<j> 2 ) + Y. I de ' K ks> i( t>a>) log (1 + Y aV (0' + i^)) 

a',s> 

-C 2 -log(l + Y 22 - , 
logY 22 (0) = -V2\m 2 \ cosh (0 - i(j) 2 ) + [ dQ' K 2,s' (0-0'- i(f> s >) log (l + Y oV (0' + ifo)) 



0-^)) -log(l + Y 32 (0-^ 



logY 32 (0) = - |m 2 | cosh (0 - ^ 2 ) + / ^'i^,' (0 - 0' - ifo) log (l + Y aV (0' + i^)) 
+ C 2 - log (l + Y 22 (0 - i~ 



Let us now examine the Y-functions appearing in the residue terms. To be specific, we will focus 
on the term involving Yi 2 (0 — in/ 4) in the first equation. For the argument of this Y- function 
we have —n/4 < Im (0 — in/ 4) < 0. As argued before, in this range of the argument no residue 
terms appear and in the multi-Regge limit the Y-function is given by 

Y 12 (0 - i^j a* e -|m 2 |cosh(e-if-i0 2 )-C 2 = g-CsgOOBh^isinM ( 52 0) 

The exponent of e 2 is positive in the chosen range of 0. Therefore, the function goes to zero in the 
multi-Regge limit and the residue containing Yi 2 gives a negligible contribution. Analogously, 
we can analyze the remaining residue terms. It turns out that they are all negligible. The same 
is actually true for n = 8 external gluons. However, starting from 9-points the residue terms 
give significant contributions in the fundamental strip, see Appendix C for explicit expressions. 
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6 Bethe ansatz 



The analysis outlined above has shown that the non-linear integral equations that control the 
n-gluon amplitudes at strong coupling simplify drastically when we take the multi-Regge limit. 
In the limiting regime we can actually neglect the integral contributions, possibly after taking 
some residue terms into account. Such a limit is well known in the theory of integrable systems. 
It corresponds to a large volume limit in which the solution of the integrable model boils down 
to solving a set of algebraic Bethe ansatz equations. 

Before we recall how the Bethe ansatz equations emerge in the multi-Regge limit, we need to 
comment a bit more on the residue terms we picked up while sending the phases <j) a to their 
limiting values ip s = (1 — s)tt/4. The appearance of such residue terms has been discussed for 
the AdSs Y-system in Appendix B of [23]. Using experience from closely related wall-crossing 
phenomena (see [36]), the authors of [23] demonstrated how residue terms can be absorbed in 
a redefinition of the Y- functions. While bringing the equations back into the standard form of 
a Y-system, 



it is necessary to introduce additional Y- functions. This implies that the new set of equations 
(6.1) might involve more than the 3ra — 15 Y- functions we started with. The complete set of 
Y- functions is enumerated by the label A,B,.... Constructing the source terms pa and the 
kernel functions Kab is part of the task one has to address while passing from a modified Y- 
system with residue terms back to a Y-system of the form (6.1). We defer a detailed discussion 
of this procedure for the AdS§ Y-system at 4> s = (1 — s)7r/4 to a future publication. Let us only 
mention that the source terms pa for the new Y-functions that are introduced while rewriting 
the equations possess the canonical form with masses tua that can be obtained from the n — 5 
masses m s of the original Y-system. 

Once we accept that the original Y-system with large phases (fi s = (1 — s)7r/4 can be brought 
into the form (6.1), we are prepared to review how Bethe ansatz equations emerge. In order to 
do so, we represent the kernel functions Kab through new objects Sab, 



As observed first by Dorey and Tateo [25, 26], upon analytic continuation of the parameters 




(6.1) 



Kab(9) 



d e \og S A b{9) 



(6.2) 



some of the solutions to the equations Ya(6) 
those solutions by an index v = 1, . . . , Na' 



1 may cross the real axis. We shall enumerate 



Y A (e^) = -l, for v = l,...,N A 



(6.3) 



When this happens, the integral on the right-hand side of equation (6.1) picks up a residue 
term since there is a pole crossing the integration contour. Hence, after analytic continuation 



the equations (6.1) take the form 



N B 

- log Y A (6) = Pa(0) + E lo § 5 ^(^ " 

S v=l 

+ E / ? A '^(^ - O lo § (i + Yfl(fi')) • (6-4) 
B 71 

If we now send the parameters of our nonlinear integral equations back into a regime where 
the integrals can be neglected, e.g. into the multi-Regge regime we explored in this note, the 
equations (6.4) become 

N B 

B u=l 

We can exponentiate this set of equations for the functions Y A (9) and insert the values 9 = o[f^ 
satisfying eq. (6.3) to obtain 

N B 

e iRk ^ = HI[S AB (0^ ~0i B) ) (6.6) 

B v=\ 

with k A (0) '■= ipA(6)/R- The parameter R and the functions k A were introduced here to help 
interpreting the equations (6.6). In our context, these equations simply determine the possible 
location of the solutions Ov to the equations (6.3). But the form of the equations coincides 
with the usual Bethe ansatz that imposes single- valuedness of wave functions for particles on a 
1-dimensional circle of circumference R. The term exp(iRk A ) accounts for the phase shift of a 
freely moving particle with momentum k A (0{/ ) when we take it once around the circle. The 
remaining factors arise from the scattering with other particles that may be distributed along 
the 1-dimensional circle. Hence, the quantities S A b introduced in eq. (6.2) are interpreted as 
a scattering matrix for excitations of some integrable system and the source terms k A describe 
the momentum. 

To make the Bethe ansatz equations (6.6) for the multi-Regge limit of the bubble ansatz more 
explicit, we need to determine the range of the label A, the source terms p A (9) and the kernel 
functions K A b(0), or rather the corresponding S-matrices. Only in the case of n = 6 external 
gluons these can be read off easily from the original Y-system. The general case will be addressed 
in a future publication. Much of the above would have remained valid if we had not brought to 
the modified Y-system back into the form (6.1). But the resulting Bethe ansatz equations (6.6) 
would have been modified as well, with the left hand side being replaced by a sum of products 
of 'phases' exp(iRk atS ). Examples can be worked out from our formulas for the multi-Regge 
limit of the Y-system with n = 9, see appendix C. Such modified Bethe ansatz equations can 
certainly be studied numerically. Nevertheless, we believe that a deeper understanding of the 
underlying integrable system requires to absorb the residue terms so that the Bethe ansatz 
equations take the standard form. 

Let us finally discuss the form of the free energy (3.12). As explained in [23] the original 
expression remains valid in the presence of large phases <f> 8 , but it needs to be rewritten in 
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terms of the Y- functions Y^. Upon analytic continuation of the parameters, the Y- functions 
can give rise to pole terms that cross the real axis. This happens precisely when the conditions 
(6.3) are satisfied. After taking the multi-Regge limit, only these pole terms survive and one 
should obtain an expression of the form 



with some energy functions = e^(0) that must be determined while passing from the modified 
Y-system with residue terms to eqs. (6.1). Hence, in order to evaluate scattering amplitudes 
of strongly coupled M = 4 SYM theory in the multi-Regge limit, we have to solve the Bethe 
ansatz equations (6.6) for the positions d\f^ of the solutions to eq. (6.3). Once these are found, 
we can easily evaluate Af rcc . 

7 Conclusions and outlook 

In this note we studied the multi-Regge limit of scattering amplitudes in strongly coupled N = 4 
SYM theory for any number n of external gluons. As reviewed above, the remainder function 
Jl( n ) [ s determined through an auxiliary ID quantum system (3.2)-(3.4) which depends on 3n— 15 
parameters m s ,(j) s and C s . The latter map in a complicated way to the 3n — 15 independent 
cross ratios that parametrize the scattering process. Our central result (5.2) identifies the 
values of the parameters in the ID quantum system that correspond to the multi-Regge limit 
of the 4-dimensional gauge theory. Since the relevant limit (5.2) involves sending all the mass 
parameters m s to infinity, the multi-Regge (high-energy) regime of the gauge theory maps to 
the large volume (low energy) limit of the ID quantum system. In such a limit, the ID quantum 
system simplifies drastically. More precisely, the non-linear integral equations (3.2)-(3.4) can 
be replaced by a much simpler set of algebraic equations (6.6). 

We find it remarkable that the computation of scattering amplitudes simplifies at both weak 
and strong coupling. In the quest for the exact S-matrix of N = 4 SYM theory, one that 
interpolates all the way from weak to strong coupling, it could therefore pay off to consider the 
multi-Regge regime as an intermediate step before addressing general kinematics. Our results 
suggests that the multi-Regge regime could be tractable even for finite coupling, at least more 
tractable than the full dependence of the remainder functions on all the 3n— 15 cross ratios. On 
the other hand, the Regge-limit imposes very strong constrains on the analytical structure of 
the remainder functions R^ n \ Explicit formulas for the Regge-limit of the remainder functions 
could therefore be an important ingredient in reconstructing the full scattering amplitude from 
more basic data. 

But before thinking about the interpolation to finite coupling, there are a few more immediate 
issues to be addressed. One is related to wall-crossing phenomena we briefly mentioned in section 
6. Recall that our limit (5.2) involves large phases 4> s in which differences \4> s — 4> s ±i\ assume 
the critical value ir/4. As we described in much detail, these large phases bring additional 
residue terms into the non-linear integral equations of the Y-system. This happens starting 
from n = 7 external gluons. The corresponding modified equations can be obtained through 
an algorithm we outlined in section 5.3. The simplest non-trivial example was spelled out 





A v=\ 
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explicitly in subsection 5.4. Following a procedure that is inspired by the study of wall-crossing 
phenomena [36, 37], it should be possible to bring the modified Y-system into the standard 
form (6.1). Our discussion in section 6 assumed that the necessary steps have been carried out 
already. But in order to determine the precise range of the labels A, B, the momenta pa(Q) an d 
the S-matrix elements Sab that enter eq. (6.6) for n > 6 external gluons, one cannot avoid a 
detailed analysis. We leave this to future research. 

Finally, we need to establish a map between the analytic continuation of the kinematic variables 
and the numbers N A of Bethe roots 9„ in the previous section. In the case of the hexagon, 
such a correspondence was determined through numerical studies. The authors of [24] found 
that the analytical continuation from the so-called physical to the mixed regime (see [24] for 
precise definitions) makes two solutions of eq. (6.3) cross the real axis. Hence, the multi-Regge 
limit of the n = 6 amplitude in the mixed regime corresponds to the energy of a doubly excited 
state in the ID quantum system at infinite volume R. This agrees nicely with the analysis in 
the weakly coupled theory. For the Regge-limit of the full 2-loop hexagon remainder function 
to be non-zero, one needs to pass into the same mixed regime that is associated with a non- 
trivial doubly excited state of the Bethe ansatz equation (6.6). It is clearly desirable to extend 
such studies beyond the case of the hexagon. Note that for larger number of external gluons 
there exist many different regimes with non-vanishing Regge-limit which probe eigenvalues of 
the BKP Hamiltonian [38] with increasing number of sites [39]. We expect that such different 
regimes are associated with solutions of eqs. (6.6) with an increasing number N A of Bethe 
roots. 
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A Explicit values of Y-functions for 7-point amplitude 

This appendix contains the complete analysis of the restrictions that the desired multi-Regge 
behavior of cross ratios imposes on the limiting values of the angles (f) s for n = 7 gluons. Some 
part of the required analysis was included and explained in section 4.3. Here we simply state 
the remaining set of formulas without further comments. Let us begin with the simplest cross 
ratios whose evaluation does not require any use of the recursion relation (3.10): 



U31 







Y, 



2.2 



-\/2m,2 cos (f>2 i q 



!>2 G 



7T 7T 

2' 2 



u 3 2 



r[l] 

2,1 



-\/2mi cos^— 4>i) 



7T 37T 

4'T 



U22 



Y. 



Ml 



2.1 



^ g-v^mi cos(^+0i} 



3tt 7T 

— 4 



- 28 - 



Since the remaining cross ratios involve values of the Y-functions outside the fundamental strip 
we must use the recursion relations (3.10). This gives: 

" -il 
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To determine the precise values of 4> s in the multi-Regge limit, we look at the ratios (2.30) 
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B Residue structure for arbitrary Y-functions 



As mentioned in the text, the residue structure for Y-functions with large shifts is intricate. In 
this section, we will demonstrate how the residue structure can be read off from the pattern 
introduced in the text. To do so, we will look at the specific example Y^ 4 '. The phases that 
appear in the Y-system equation are (f>i = and (p2 = — vr/4. Since we need to evaluate our 
Y-function at 9 = —in, the quantity 9 — i(pi crosses the poles —ij, — i\ while 9 — icj)2 
crosses —ij, — if • The resulting residue structure therefore reads 
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Figure 5. Residue structure of Y"J j 4 \ For simplicity, the residue dependencies of the other Y-functions 
arc not shown. 



Graphically, this residue structure can be represented as in figure 5, which remains valid for 
yj j 4 ' . Of course, some of the Y-functions that appear as residues themselves receive corrections 

[—31 

from residue terms, for example Kj :2 ^ n ^ ne aDOVe example. However, it should be clear from 
the graphical representation that the residue term with the highest s' that can contribute to 



a Y-function Ya}s with k negative is given by the intersection of the diagonal Y^^+j with the 
diagonal Y^ • . 



[k+i] 



C 9-point multi-Regge limit 

Here we present the equations governing the Y-system in the multi-Regge limit. In the lower 
half of the fundamental strip Im# 6 (— 7r/4, 0) the equations are actually obtained simply by 
dropping the integral contributions from the original expressions (3.2)-(3.4). For the upper half 
of the fundamental strip Im9 G (0,7r/4), some residue terms survive in the multi-Regge limit. 
The relevant equations take the form log Y ajS = p a)S with the right hand side given by 

P2±i, s (d) = -m s cosh (9 - i<f) a ) ± C s 

P2,s(8) = -V2m s cosh (9 - icj) s ) 

for the first two values s = 1, 2 of the parameter s. Recall that the corresponding phases cj) s are 
given by (j>± = and (p2 = — vr/4. All the remaining functions contain residue terms. These are: 

P2±i,3 (9) = -m 3 cosh (e + 4) ± C 3 + log ( 1 + e ±C4- m 4cosh(tf+ 4f ) 



P2,3 (0) = -V2m 3 cosh (6 + i|) + log (l + e ^»Mcosh(0+ lf ; 



P2±i,4 (^) = — w-4 cosh [ 9 + i 



.3vr 



± C 4 - log 1 + e 



-V2m 4 cosh(6»+«-|) 



P2A (0) = -V2m A cosh (0 + - log ((l + ^-m^e+i^ ^ + e _ C 4-m 4 cosh( e+if ) 

As mentioned in the text, these equations could be used as the starting point to find the solutions 
of equation (6.6) numerically. If we use the above functions p a ,s-, however, the left-hand side of 
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eqs. (6.6) involve sums of products of exponentials such as 

e J>2,3(0) = e -v / 2m 3 cosh(6»+jf ) A + e - v / 2m 4 cosh(0+i§)^ 

In order for the Bethe ansatz equations to take a more standard form one needs to work with 
a larger set of momenta pa and the corresponding S matrices as discussed in the main text. 
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